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Abstract 

Using an energy-momentum tensor for spinning particles due to Dixon and 
Bailey-Israel, we develop the post-Newtonian approximation for N spinning 
particles in a self-contained manner. The equations of motion are derived 
directly from this energy-momentum tensor. Following the formalism of Ep- 
stein and Wagoner, we also obtain the waveform and the luminosity of the 
gravitational wave generated by these particles. 
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I. INTRODUCTION 



The possibility of detecting gravitational waves from coalescences of compact binaries in 
the near future by LIGO [|TJ and VIRGO laser interferometric detectors has aroused quite 
a lot of research in this area. One of the main tools in analyzing these situations and in 
obtaining the gravitational waveforms generated is the post-Newtonian approximation ||. 
As one gets closer to the coalescence, higher order terms in this approximation have to be 
taken into account to get results accurate enough for comparison with future observational 
data. 

When higher order terms are considered, spin and higher multipole effects may be of 
significance, especially for neutron stars and Kerr black holes. For these compact objects, it 
is estimated that the spin-orbit and spin-spin effects are of (post) 3//2 -Newtonian and (post) 2 - 
Newtonian order, respectively ||. Therefore, to get accurate results to these orders, spin 
effects must be included. 

A number of authors have studied these spin effects in ||-||. In 0, both the spin-orbit 
and spin-spin effects are considered in details using the multipole formalism of Blanchet, 
Damour, and Iyer (BDI) J?|||. The BDI formalism is rigorous, but rather complicated in 
which each body is regarded as a spherically symmetric, rigidly rotating perfect fluid. 

In this paper we would like to introduce another approach which is simpler and self- 
contained in the sense that the equations of motion are also derived along the way. In this 
approach we still regard the bodies as point particles, but with structures such as spin. 
The central quantity here is the energy-momentum tensor for these particles. In general 
spacetimes this energy-momentum tensor was first introduced by Dixon |J, and was later 
elaborated on by Bailey and Israel (BI) |l(J . In [|Kj , the Lagrangian density of a particle is 
expanded in a covariant Taylor expansion about another world-line, presumably the center- 
of-mass of some extended object. Higher order terms in this expansion correspond to spin, 
tidal effects and so on. From this Lagrangian one can obtain the energy-momentum tensor 
in the usual fashion. Therefore, in this formalism one can include structures such as spin, 
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quadrupole moment and so on in a systematic manner. In this paper, however, we shall 
concentrate only on the spin effects. 

With this energy-momentum tensor for particles with spin, one can derive the equations 
of motion for spinning particles, that is, the Papapetrou equations, by requiring this energy- 
momentum tensor to be symmetric and conserving. This will be done in the next section. 
In Section III, we develop the post-Newtonian approximation for N spinning particles using 



this tensor. We follow closely the procedures in ||11|| . As an illustration of how our approach 
works, we choose the simpler case that the spin-orbit and spin-spin effects are both of post- 
Newtonian order. In this case, we are not dealing with compact objects, but ordinary ones. 
The equations of motion for positions and spins are expanded up to (post) 2 -Newtonian 
order. In Section IV we use the formalism of Epstein and Wagoner (EW) |~2| to obtain the 
waveform and the luminosity of the gravitational wave generated by the motion of these N 
spinning particles up to post-Newtonian order. This can be considered as an extension of 
the result of Wagoner and Will Jl3| to spinning particles. Conclusions and discussions are 
given in Section V. 



II. ENERGY-MOMENTUM TENSOR 

In ||10|| , BI devised a general method to derive the energy- momentum tensor T^ v for 
particles with structures such as spin. First the Lagrangian for a charged point particle 
is expanded covariantly about another world-line. For the particle to represent part of an 
extended object, this world-line is usually chosen to be the center-of-mass of the object. 
However, their formalism is rather general that one can actually choose any convenient 
world-line. The higher order terms in this expansion can be identified with the structures of 
the particle corresponding to higher gravitational and electromagnetic multipole moments. 
With this Lagrangian expansion, one can obtain the equations of motion and T^ u by the 
usual variations. 

Here we concentrate on particles with the extra structure of spin, which is related to the 
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gravitational dipole moment. In this case, T pu for a spinning particle is given by (10 



(2.1) 



where 



rpflU 

(can) 



dr^v u 5\x- z(t)), 



gv-vp — (S pup + S ppu + S pup ) 



guup 
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drS pu v p 5 4 (x-z(r)), 



(2.2) 

(2.3) 
(2.4) 



and p p is the canonical momentum, the spin tensor. T p can ^ is therefore the canonical 
energy-momentum tensor, and B pup is the Belinfante tensor due to spin. The equations of 
motion can be derived in the following manner. First we require T pv to be symmetric in /i 
and u, T pu = T up , which gives 



Using the formulae, 



dr (p p v u - p u v p )5 4 (x - z{t)) - V„ 



drS pv v p 8 4 (x - z(t)) 



0. (2.5) 



duV^ = V=g~ T p 



pal 



(2.6) 



and 



j p d,5\x - z{t)) 



-^(x-z(r)), 



the derivative becomes 
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drS pu v p 5 4 (x - z{t)) 



DS pu 
dr——5 4 (x - z(t)), 
Dt 



where we have ignored surface terms. Here 

DS pv 



Dt 



= v a V a S pv 

= — + v a r p ; p s pl/ + v a r» ap sw. 



(2.7) 



(2.8) 



(2.9) 



Therefore, Eq. (|2.5|) becomes 



1 



DS» V 
Dt 



dr 



DS^ 
Dt 



5 a {x-z{t)) =0 



(2.10) 



This is the first set of the Papapetrou equations, the equations of motion for spinning 
particles. 

Next we require that T^ u is conserving, V M T M!/ = 0, which gives 



(2.11) 



where we have used the symmetry property of T^ v . Using the formulae Eqs.( |2.6|) and 
again, we have 

1 



XI T vti — V 



-g 



drp u v fl 5 4 (x -z(t)) 



(2.12) 



For V „V ' P B UW ', we use the identity 



(2.13) 



for the tensor a/3 "' and the symmetry properties of the Riemann tensor Rnua/3 to get 



(2.14) 



Finally, Eq. fl2.11 ) becomes 

1 



-g 



dr 



Dp u \ 1 



5\x-z(t)) = 



Dp" 1 

— — = — R v w <tS P(T v>*. 
Dt 2 w 



(2.15) 



This is the second set of the Papapetrou equations. If the spin tensor vanishes, 

Dp v 



Di 



0. 



(2.16) 



which is just the geodesic equation for a point particle. Also from Eq.( ^.10[ ), we see that the 
canonical momentum can be written as 



p» = mv^ -v v ——. (2.17) 
Dt 

The spin tensor is anti-symmetric in /i and u, giving six independent components. 
Three of them can be eliminated by the "spin supplementary condition" , 

S^v v = 0. (2.18) 

For fi = i, 

S i0 v + S ij Vi = 



S" 1 = M- S lJ . (2.19) 
\v J 

The fi = equation is then satisfied automatically. The remaining three independent 
components of can be written vector, 

S i ^ I e *i^ifc ; (2.20) 
which is the spin vector of the particle. 



III. POST-NEWTONIAN APPROXIMATION 

In this section we develop the post-Newtonian expansion for N spinning particles using 
the energy-momentum tensor T^ v introduced in the last section. We follow closely the 



procedure in (TT| . In the post-Newtonian approximation, it is assumed that the velocities of 



the particles are small, 

w~e 1/2 <l. (3.1) 

Moreover, the gravitational potential <fi is comparable to the kinetic energy, 

~ v 2 ~ e. (3.2) 

Here we take the spin S l to be of order e 1//2 , the same as velocity. Under this assumption, the 
spin-orbit and spin-spin effects are both of order e. This is valid for ordinary situations but 



not compact objects for which these effects will come in only at higher orders. In this paper 

we wish to treat the former case, which is simpler, as an illustration of how our approach 
works. 

To begin we expand the spacetime metric, 

<?oo = -l + <?£ ) + <?$ + •••, (3.3) 

9(H = 9 { k ) + ---, (3-4) 

9ij = &ij + 9if + ■■■, (3-5) 



where 



9$ = -20, (3.6) 
gS = -2(^ + ^), (3.7) 

g& ] = 0, (3.8) 

gff = -2<t>8ij. (3.9) 

Here the superscript denotes the order in the post-Newtonian expansion, for example, 
is the (post) l//2 -Newtonian order of #oo- Through the Einstein equations the potentials 0, Q, 
and ifi can be expressed in terms of the various components of T^ v in the appropriate order, 

<p(x,t) = - 1 rfv (T °y /; ff (o) , (3.io) 

J \x — x \ 

Q(x,t) = -A /^V (T °y /; ^ )(1) , (3.11) 
J \x — x'\ 

1,(2, t) = ^ - ! d 3 x> [(T 00 (f ,t))( 2 ) + (T\x',t))^\ , (3.12) 



dt 2 

where V 2 x = 4>- Note that we have used the harmonic gauge [11 



gn% = 0. (3.13) 

in deriving the above expressions. 

As for T^ v , we use the one given in the last section for N spinning particles, 

TflU = T {can) + VpS""', (3.14) 



where 

Tfel) = -^EpW 3 ^-^), (3.15) 

B^p = --(S^p + S ppu + S pup ), (3.16) 
2 

S w = _L = <£s!?v> A 6 3 (a!-2 A ), (3.17) 

V # A 

= m A v A - v A u ^ A , (3.18) 
Dt a 

and A = 1, N. Assuming that the spin part is of the order of e 1 / 2 , we can expand T 00 , 
T 0i , T i0 and T ij , 

(T 00 Y 0) = J2 m AS 3 (x-x A ), (3.19) 

A 

(T oo )(2) = ^ mA ( 0A + 1 2 )5 3 (f _ fA) _ ^ e ^M^ 3 (f - X A ), (3.20) 

A A 
(T 0i)(l) = £ m ^ 3 (£ _ f a) + 1 £ e Ufe^ 9 . 5 3 (f _ f a) 
A A 

= (r°) (1) , (3.2i) 

(T ^ = E ™a(0a + ^iK5 3 (f - x a) + E e*M (^j 5 3 (f - x A ) 

+ E e ^ (^) ( } ^ 3 (5 - 2a) + E ^0a^ 3 (* - X A ) 
-\ E ^>^5 3 (:r - £ A ) - i E ^> A ^ 3 (f - £ A ) 

A A 

/ jqk \ ( 2 ) 

= (r°)( 3 ) + E^M (^J 5 3 (^-^a), (3.22) 
(t^)( 2 ) = E^> A 5 3 (f - x A ) - ^E^' fc (^J £0 

+ ^ {lkl v°is A d k 5\x-x A ), (3.23) 

A 

(r*) (4) = E ^(0a + ^»^ 3 (* - ^a) - ^ E (^J * 3 (* - £0 

+ E t {lkl <s A ^ 5 3 (f - f A ) + E ^ ^ 5>^ 3 (f - x A ) 



+ E e -(#) (2) ^ 3 (-^) + E^^)^-^ 
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/ jql \ (2) / JQfc \ ( 2 ) 

+ E ^ (^J ^ 3 (* - *a) - £ (^J <M 3 (* - *a) 
+2 £ e^vHs^AdkS 3 ^ - x A ) - E e (ifem ^>5^ 3 (£ - £0 

+E^M(#M* s (2-aU) 

+2 £ e iU v k A S l A {d?<l> A )P(2 -x a )-2J2 ei kl v A S A (d£cf> A )5 3 (x - x A ) 

A A 

~\ E ^fiiC^U)* 5 ^ -x A ) + W e^ l S l A (diU^\x - x A ), (3.24) 

A A 

where we have used the convention for symmetrization of indices, 

A (i B j) = -(A i B j + A j B l ) (3.25) 
2 



From this expansion of T"", the potentials (j), Q, and V in Eqs.(glQ|), (pTD and ( gig ), 
respectively, can be evaluated to be, 

= -E^% (3-26) 



.4 



^A^A . n {x - X A ) X SU 



Yf-^a| x |#-aul 
^ _ m A m B 1 m A (x - x A ) m B (x A - x B ) 

A B^A ft - \XA ~ Xb\ 2 AEh£A \$ ~ %a\ \xa~X B \ 3 

| 1 m A [v A ■ (x -x A )f 2 y m A v\ 2 y (x - x A ) ■ (v A x S A ) ^ ^ 

2* ^ I rft <-y> , I 3 ■ ^ I nn y , I ■ ^ I rv* np , I 3 

A \x xa\ a I A A I A | 

The equations of motion for xa and Sa can be derived, as discussed in Section II, by 
requiring T^ u to be symmetric and conserving. For T^ u to be symmetric, we have the 
conditions, 

T ij = T ji_ (3 3^ 

On the other hand, for T^ v to be conserving, we have 
V^ = 

=>> |. t qo + + 2 r0 T 00 + 3T° 0l T i0 + Ti T 00 + T° H T ij + T^T = 0, (3.31) 



and 



V M T^ = 




dt T 0i + djT ij + T° Q0 T 0i + r° oj T ij + T k kQ T° l 
+T k kj T» + P nn T 00 + 2TlT 0j + T] h T jk = 



oo 



jk- 1 



[3.32) 



The lowest order of Eq. (|3.29|) is satisfied automatically because (T 01 )^ = (T* )^ as given 
in Eq.flOlD. From Eq. flOOp , we have 



(rpij\(2) __ (rpji\{2) 
' - \ (2) 

<ISa 



dt 



0. 



(3.33) 



which is the equation of motion for spins to the lowest order. From Eq. ( P-31|) , we have 



' d_ 
dt 



\ (o) 

T 00 + <9,(T° 4 ) (1) =0, 



(3.34) 



which is also satisfied automatically, and so there is no new constraint. From Eq.( |3.32| 

(i) 



9_ T oi\ + 9i(r ii)(2) + (r 00 )( 2 )(T 00 )(°) = 



dVA 

dt 



(2) 



-E 



m B r A B 



BjtA r AB 

where we have denoted xa — %b by tab- This equation is just the Newtonian equation of 
motion for point particles. The various components of the Christoffel symbol (T^py' can 



3 ' 



(3.35) 



be found in 

Now we go to the next order. From Eq. ( |3.29| ), we have 



(T 



0^(3) = (j-i0\(3) 



(3.36) 



which is satisfied automatically as seen from Eqs. (|3.22j ) and (|3.33| ). Then from Eq.( [3.30| ), 
we have 



(T 



m(4) 



m(4) 
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dS* 
~dt 



(i) 



^2 (^fj + va • V A <pjij S A - {S A ■ v a )Va4>a 

+2(S a ■ Va<Pa)v a - \{S A ■ V a )Ca + ~Va(& ■ CO 
m Br A B ■ (v A ~ 2v B ) g ^ m B [(v A - 2v B ) ■ S A ) 

1^ 13 A 2^ a r AB 



B^A 



AB 



B^A 



AB 



, m B (jAB ■ Sa) ,-, — * \ Sa X S B 



B^A 



AB 



B+A ' AB 



t Q s-^ r AB ■ (S A x S B ) ^ , o (tab -Sa),-, xx ^ v 
+3 2^ rAB + 3 5 ( r AB x S B ), 



Bj^A ' AB Bj LA 1 AB 

where we have used the results in Eqs.( |3.26| ) to (|3.28| ) for the potentials. 
Next from Eq.(p3T|), 



^ rriOO 

dt 



(2) 



+ ^( T 0^(3) + 2(r O o) (3) (T 00 ){ 0) 
+3(r 4) (2) (T iO )( l) + (r y(3) (T 00 )( 0) + (r , )(2 ) (T o,. )( l) = Q ; 



(3.37) 



(3.38) 



which is satisfied automatically because of the Newtonian equation of motion as in Eq.( [3.35| ). 
Then from Eq. flO^) , 

\ (4) 



d_ 

dt 



T 0i + dj(T 1 ^ + (r° )( 3 )(T 04 ) (1) + (r^) (2) (^') (2) 



+ (r * o) (3) (T 0* )( l) + (r y(2) (r , ){ 2) + (no) (2) (T 00 )( 2) 

+(rj )0 )( 4 )(T 00 )(°) + 2(rL) (3) ( r ° i ) (1) + ( r i*) (2) ( ri *) (2) = o 



(3.39) 



which gives 

/ dv A \ 
{dt 



("1) 



3v A 



dt 



v 2 a (Va<Pa) + 4:Va(va ■ Va0a) 



-40 A (V ' A (Pa) - VaVu - -qT- + v a x (Va x (a) 



m A \ 



dt 



S A x V A (v A • V a <Pa) 



m 1 



(v A x S A ) ■ Va</>a 



+ 



2rriA 



S A ■ (V A x (a) 



E 



m B r AB 



B^A r AB 



4£^+ £ 



C^A T AC C ^A,B T BC 



BC 
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+5 v A + Av A ■ v B - 2v B + 



tab 



2 .2 



r AB 



7 ^ ^ mcf B cm B , v-^ m B r AB • (4v A - 3v B ) „ , 
L ~3 — + E ~3 (va-vb) 

* B^AC^A,B ' BC ' AB Bj l A ' AB 



-6E 



tab 



[tab x (v A - v B )\ ■ {m A S B + m B S A ) S A ■ S B 



Bj l A m A r AB 
5 (tab ■ S A )(r AB ■ S B ) 



AB 



■ AB 



+ 



AB 



1 



-3 E — V - [(Tab ■ S A )S B + (f AB ■ S B )S A 

B^A 

m A r AB 

- E — ~3~ (va - vb) x (Am A S B + 3m B S A 

B ^A m A r AB 



+3E 



1 



-TAB ■ (VA - V B 



tab x (2m A S B + m B S A ) 



(3.40) 



B -£ A m A r AB 

Hence, Eqs. (|3.33|) , (|3.35|) , ( |3.37|) , and ( |3.40| ) form a complete set of equations of motion 
for x A (t) and S A (t) for spinning particles to the (post) 2 -Newtonian order. This can be 



considered as an extension to the well-known Einstein- Infeld-Hoffmann equations [13] for 
ordinary particles. 

Using these equations of motion, the expressions for the energy-momentum tensor T^ v 
in Eqs.(|3.l9|) to (3.24) can be simplified to 



(y00^(0) 
(j-00\(2) 

/J-OiN(l) 



A 

T, m A(4>A + \v A )5\x-x A ) -Y,£ ijk v A S k A di5*{x-x A ) 

A A 



(T 



iO\(l) 



1 



E m A v A S 3 (x - x A ) + - E e^S^S^x - x A ) 

A 1 A 



(3.41) 
(3.42) 

(3.43) 



(T 



on (3) 



(T 



i0\(3) 



1 



E ^a(0a + ^i)^ 3 (^ - £4) 



+ e l > k S A (d£<p A )5 3 (x - x A ) - E e ijk S 3 A cl> A d k 5 3 (x - x A ) 



-^ kl ^ k A S l A d j 5\x-x A ) 



-^ kl v\v k A S l A d 3 8\x - xa) 

Z A 



(T 



1 -J * ■) + E ^vH&dkPia! - x A ) 



J2m A v A v A S"{:r - .r A 
A 



(3.44) 
(3.45) 
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Cn (4) = £ m A ^ A + \v A )v\v A b\x - x A ) 

A / 



+2 E e^^0 A (9 fc 5 3 (f - f A ) - 4 £ e^^A^A^x - x A ) 

A A 

+ E e im v k A S A (&>2<pA)5 3 (x - fx) - E e (ifcm ^S>^<5 3 (f - x A ) 

-~E ^a(^U)5 3 (x - f x ) + i E ^(#CiV(a? - *a), (3-46) 

A A 

and the expressions for the potentials 0^ and Ca are those given in Eqs. fl3.26]) and fl3.27| ). 
These components of the energy-momentum tensor T^ u to various post-Newtonian order 
will be used in the next section to find the gravitational wave generated by a system of N 
spinning particles in the post-Newtonian approximation. 



IV. GRAVITATIONAL WAVES GENERATION 



In this section we follow the works of EW [|T2| and Wagoner and Will ||T3| to consider the 
gravitational wave generation of N spinning particles in the post-Newtonian approximation. 



As in ||12|| , we derive the waveform and the luminosity of the gravitational waves up to 
post-Newtonian order. 

The waveform of the gravitational wave in the radiation zone is given by fl2 l 



2_ a_ 2 

Rdi 2 



h TT = ^^j I ij (t - R) + n k I^ k (t -R) + n k nd iik \t - R) 



TT 



where TT denotes the transverse-traceless part, n k = Rk/R, and 



1*3 = I T 00 xVd 3 X, 

jijk = J (2 T ^ i x j) x k - r 0k x l x j )d 3 x, 



Jijkl = / T H x k x l d 3 x 



where t^ u is the total energy-momentum tensor, 



which includes the matter part T^ v and the gravitational part t^ v ' . 



(4.1) 



(4.2) 
(4.3) 
(4.4) 



(4.5) 
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Here, for spinning particles, we use the energy-momentum tensor T^ v discussed in the 



last section. While t^, up to post-Newtonian order, is given in ||13||. Therefore, 



r 00 = £ ™a(1 + <Pa + \v 2 A )5\x - x A ) - E e i]k v\S k A d t 5\x - x A ) 

A 1 A 

-J-(4^ + 39i^), 
r° l = J2 m A v A 6 3 (x -x A )-W e^ k S A d k 5 3 (x - x A ), 

A 1 A 



t 13 = m A v A v A 5 3 (x - x A 



) + Y^ m S l A v i A ) d k 5 3 (x-x A ) 



(4.6) 
(4.7) 



-2d i( j)dj<j) - 44>did j( f) + %(40<9 2 + 3d, 



'k<pv k i 



(4.8) 



up to post-Newtonian order, where the last terms in Eq. (|4.6| ) and Eq. (|4.8| ) come from the 
contributions of t 00 and t lJ , respectively. Then the integrals in Eqs. (|4.2| ) to ([4.4]) can be 
evaluated, 



/^^m^4 1 + 



1 



1 



: V A ~ 



E 

B^A 



m B 



2 * 2 — r AB 



I l3k = E rn A (2vMx k A ~ v k A x A x\) + 2 £ e k ^4s l A , 



(4.9) 
(4.10) 



1 ^ ^ m A m B r l AB r AB 



Y, (e {ikm v A ] x l A + e (i 



A B^A 
lm„J)„k 



rAB 



1 



S ki + (-r AB r AB + 6x A x l A ) 

r AB 



i^JIrfk \ cm 

v A x A \ o A 



(4.11) 



In obtaining these expressions, we have discarded terms with vanishing transverse-traceless 
parts, of which some are actually divergent, because those terms will not contribute to the 
generation of gravitational waves Hl3| . For example, 



tJ/TT 



0t) 



TT 



0. 



0. 



(4.12) 
(4.13) 



for some function fK Also we have used the formula, 



%dj(f)x k x l d 3 x 



EE 

A B^A 



Tcm A m B r AB r AB 
3r AB 



^ H + ~j2~(2r AB r l AB — Qr AB x A + Qx A x L A ) 



( k J) 



AB 



(4.14) 
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which can be proved [|T^] by ignoring terms which have no transverse-traceless parts in a 
similar way. 

Putting these together the waveform of the gravitational wave, up to post-Newtonian 
order, is 



2 d 2 J _ i , / _ 1 2 1 _ m B 

k TT = 1^^AX A X\ \ 1 - U ■ V A + -V A - - ^ -jg 



Rdt 2 



A \ " B+A 

+2^m#S4(fi • xa) + ^Zm A v l A v A (h ■ x A f 

A A 



1 ^ ^ m A m B r AB r AB 



-—y y 

12 ^ ^ 



A B^A 

.(*/ 



tab 



1 + -=— (-(n ■ r AB f + 6(n • x^) 2 



AS 



+2$>2(ijU x S A f-2j2x { l(n x ^ - 2 ^(n • x S A )» . 

AAA ) XT 

(4.15) 



Hence, together with the equations of motion in Eqs.( |3.33"D , ( |3.35| ), ( p.37| ) and ( |3.40| ), one 
can obtain the gravitational waves generated by N spinning particles. 
The total luminosity of the gravitational wave is also given in [|I^| as, 



-QAfijkM ikj + llMijM ijkk - 2AMijHikkj) + •••), 



(4.16) 



where the bracket denotes an average over several characteristic wavelengths, and 



Mi 



^3 I 1 ijkx-k rl 



Up to post-Newtonian order, L can be expanded to 



ij ikkj J 



(4.17) 



(4.18) 



Using the results for I % \ P jfc , and I^ kl in Eqs.( |4.9|) , (|4.10 ), and ( 4.11 ), respectively, one can 
obtain the various orders of Nijki-km as follows. 

1 



AT (0) = — 
13 dt 3 



J2 m A \x\x\- ^x Aj ( . 



(4.19) 
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(4.22) 

This completes our consideration of the gravitational wave generated by the motion of N 
spinning particles to the post-Newtonian order. 

V. CONCLUSIONS AND DISCUSSIONS 

In this work we consider a method to develop the post-Newtonian expansion for both the 
equations of motion and the generation of gravitational wave of N spinning particles. In this 
approach the energy-momentum tensor T^ v for spinning particles, due to Dixon and BI, is 
introduced first. Using this T^ v , one can derive the Papapetrou equations by requiring it to 
be symmetric and conserving. By the formalism of EW, we can then obtain the waveform 
and the luminosity of the gravitational wave generated by the motion of these spinning 
particles. 

This approach is straight forward and much simpler than the BDI multipole formalism 
since we still regard the particles as point-like. It is self-contained in the sense that once we 
write down T^, the equations of motion can also be derived from it. Another merit is that 
the BI formalism is very general so that structures other than spin, such as tidal effects, 
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can also be taken into account systematically by introducing terms corresponding to these 
structures in T^ v . 

However, there are several points in our calculation that we must be careful with. First, 
the EW formalism is less rigorous than the BDI formalism in which integrals like those in 
Eqs. (|4.2|) to ( |4.4j ) have divergent parts. Fortunately, these divergent parts have no transverse- 
traceless parts so they will not contribute to the generation of gravitational waves. We 
therefore expect the two formalisms to give the same results. 

Another point to note is that the harmonic coordinate condition used in Section III, 



as well as in [fLljj , is not the same as the gauge condition used in the EW formalism. It 



is nevertheless shown in [15| that, at least to the post-Newtonian order that we are con- 
sidering here, the extra terms coming from these different gauge choices all have vanishing 
transverse-traceless parts. Therefore, there is no inconsistency in our investigation in this pa- 
per concerning gauge conditions. On the other hand, if we want to extend our consideration 
to higher post-Newtonian orders, we must be more careful with gauge choices. 

In our calculation we have assumed that spin S is of post-Newtonian order, same as 
velocity v. Then the spin-orbit and spin-spin effects come in both at the post-Newtonian 
order. Since these effects for compact objects come in at even higher orders, we are really 
dealing with ordinary objects. We have chosen this simpler case just as an illustration of how 
our formalism works for spinning particles while staying at the post-Newtonian level. Of 
course, the compact case is more interesting because of the possibility of detection of their 
coalescences by LIGO and VIRGO in the near future. We plan to consider that situation in 
a separate work. 
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